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Abstract 

The Liquid Xenon Gamma-Ray Imaging Telescope (LXeGRIT) is the first re- 
aHzation of a liquid xenon time projection chamber for Compton imaging of MeV 
7-ray sources in astrophysics. By measuring the energy deposit and the three spatial 
coordinates of individual 7-ray scattering points, the location of the source in the 
sky is inferred with Compton kinematics reconstruction. The angular resolution is 
determined by the detector's energy and spatial resolutions, as well as by the sep- 
aration in space between the first and second scattering. The imaging response of 
LXeGRIT was established with 7-rays from radioactive sources, during calibration 
and integration at the Columbia Astrophysics Laboratory, prior to the 2000 balloon 
flight mission. In this paper we describe in detail the various steps involved in imag- 
ing sources with LXeGRIT and present experimental results on angular resolution 
and other parameters which characterize its performance as a Compton telescope. 



Introduction 



The Liquid Xenon Gamma Ray Imaging Telescope (LXeGRIT) is a prototype 
of Compton telescope (CT) based on a liquid xenon time projection chamber 
(LXeTPC), which combines good energy resolution with imaging of individual 
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MeV 7-ray interactions with submillimiter position resolution (Sec. [T]). A CT 
images 7-ray sources in the energy range from few 100 keV to more than 
10 MeV, reconstructing the direction of individual 7-rays through Compton 
kinematics. The scatter angle on a free electron {(p) is given by the Compton 
formula 
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where mgC^ = 0.511 MeV, Eq is the initial energy and Ei is the energy of the 
scattered 7-ray. Fig. [T] helps illustrating the principle of Compton imaging. A 




Fig. 1. Principle of a Compton telescope 

7-ray from a source at position (xq, ^0) Compton scatters in the plane Di by 
the (true) angle ip and is stopped in the plane D2; is estimated from Eq. [TJ 
The time sequence of the interactions needs also be known (Sec. Interaction 
positions are measured in both planes and give the direction of the scattered 
photon (x, Ip) The direction of an individual 7-ray is then determined to an 
event circle with radius around the direction (x, V")- This ambiguity stems 
from the non-measurement of the direction of the scattered electron. After col- 
lection of many source events, intersection of all event circles defines the source 
position (Fig. [2]). While back- projection of event circles on the sky is useful in 
visualizing the basic measurement principle, it is not the optimum method to 
perform Compton imaging. In fact, it doesn't include the probability distribu- 
tion for scatter angles da /dip for photons from a source at given location and 
energy, i.e. the frequency of occurrence of event circles with a given opening 
angle. Fig. [2] illustrates how this additional information can be exploited in 
a 3D data space consisting of (x, ^z^,*^), where a point source defines a cone- 
like structure with half opening angle of 45° centered on the source position, 
with a density along the p dimension given by the Klein-Nishina cross section 
convolved with factors that result from the detector geometry and detector 
thresholds. Compton imaging of 7-sources is discussed in detail in Sec. [3] and 
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Fig. 2. Event circles for a point source (left) and the corresponding response function 
in the 3D imaging data space of an ideal classical Compton telescope (right). A point 
source describes a cone with half-opening angle of 45°. The density along the scatter 
angle dimension corresponds to the Klein-Nishina cross section and describes the 
density of circles of radius in the event circle representation [l] . 

examples are given in Sec. IH 



1 The LXeGRIT Compton Telescope 



For a detailed description of LXeGRIT and its performance as an imaging 
calorimeter, see [2|3lH] . In this section we summarize those aspects most rele- 
vant to Compton imaging. 

The fiducial volume of the LXeTPC is a box with dimensions 18.6 x 18.6 x 7 cm^ 
(in the x, z coordinates respectively), filled with high purity LXe. At a 
temperature of T = —95° C the density of LXe is 2.85 g cm~^ [5] and the 
attenuation length for 1 MeV 7-rays in LXe is ~6.2 cm [6]. When a 7-ray in- 
teracts in the fiducial volume, both scintillation light and ionization charge are 
produced efficiently, with W- values of VTph = 24 eV and = 15.6 eV [S]. 
The VUV (178 nm) scintillation photons are detected by four photomultiplier 
tubes (PMTs), which provide the event trigger and the initial time, to. The 
PMTs Lil are coupled to the LXe volume via quartz windows. The ionization 
electrons drift under an applied electric field of 1 kV / cm, inducing a signal on 
two parallel wire planes, after passing a Frisch grid. There are 62 wires in each 
plane and the pitch of the wires is 3 mm. The location of the hit wire(s) in the 
two wire planes provide the x and y coordinates in the TPC reference frame, 
while the time, measured starting from to, gives the interaction depth {z co- 
ordinate). The wires are transparent to the drifting charge, which is finally 
collected by one of four independent anodes, and the amplitude measures the 
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Fig. 3. Schematic of the LXeTPC readout structure (not to scale) with correspond- 
ing Ught trigger and charge pulse shapes. From [3]. 

energy deposited in the interaction. A schematic of the readout structure and 
light trigger of LXeGRIT is shown in Fig. [3l 

A reliable Monte Carlo (MC) simulation of the detector has been developed 
and tested to produce accurate results up to several MeV[3j. For a generic 7- 
ray event, the LXeGRIT output is (Ei, xi, yi, zi); (E„, x„, y„, z„) where n 
is the event multiplicity, i.e. the number of detected interactions once the finite 
spatial resolution and energy threshold (~150 keV for a single interaction) are 
accounted for. 



1.1 Position resolution 

The interaction location in the LXeTPC is obtained from the wire signa- 
ture. With a wire pitch of 3 mm, the spatial resolution is 0.87 mm (rms = 
3 mm/v^T2) in a; — if only a single wire signal is detected for each coor- 
dinate. For LXeGRIT this is the typical case for energy deposits <0.3 MeV 
[9]. The reconstructed image of a collimated beam of 7-rays (0.662 MeV) pho- 
toabsorbed in the sensitive volume is shown in Fig. |H both in the x — y and in 
the X — z views. The ^^^Cs source was located above the TPC, collimated to a 
beam with 3 mm diameter in the x — y plane by a lead collimator 15 cm thick. 
The position resolution in the z-coordinate is estimated to be about 0.25 mm 
(1 sigma). It is derived as the difference in drift times measured independently 
on the X and y wires for the same interaction, as shown in Fig. [5]-/e/t. As- 
suming two independent measurement with Gaussian errors (neglecting the 
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uncertainty in the common trigger), a = 0.35mm/-\/2 ^ 0.25 mm. The z- 
distribution of photoabsorbed events from the same colhmated ^^''Cs source is 
shown in Fig. fright, well reproducing the expected exponential attenuation. 
This spatial resolution fulfills the requirement of a fine grained CT, where a 
typical separation between interactions is in the few cm range and the linear 
extension of the charge cloud due to MeV energy deposits is typically less than 
1 mm. 




Fig. 4. Image of the 7-beam from a collimated Cs source on top of the TFC. Left: 
projection in the x — y plane. The coordinates have been re-defined in order to have 
the source image centered at x = 0, y = 0. Right: side view {x — z plane) of the 
same 7-beam. 




Fig. 5. Left: difference between the z-position as determined from the x and y wires. 
Right: z-distribution for a colhmated ^'^^Cs source sitting on top of the LXeTPC. 
Events in the full energy peak have been selected; superimposed {dashed line) the ex- 
ponential attenuation for 0.662 MeV photons. The slight discrepancy at 2; ~ 60 mm 
is mainly due to the z dependent light trigger efficiency, which is not corrected for. 



5 



1.2 Energy resolution 



The energy response is determined to be linear over the energy range from 
0.5 to 4.4 MeV, covered by the available calibration sources 0, and is shown 
in Fig. O-Ze/t. The energy dependence of the energy resolution is shown in 
Fig. fright and is described as 



AE/E{FWHM) 




where Pi accounts for the intrinsic energy resolution from the statistics of 
charge carriers in LXe; 1/\/E reproduces the energy dependence expected 
from Poisson statistic; P2 accounts for contributions which are independent of 
energy [noise term). The noise term parameterizes electronic noise, errors in 
fitting the anode wave function, shielding inefficiencies of the wires structure 
etc. In practice it is well described considering the electronic noise only (63 keV 
FWHM, which corresponds to about 1000 equivalent noise charge, as in [4|). 
The intrinsic energy resolution is 8.3% at 1 MeV, in good agreement with 
previous measurements in LXe at the same drift field (1 kV/cm), from gridded 
ionization chambers of much smaller fiducial volume |10j . 
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Fig. 6. Left: linearity plot for ADC channel vs. energy in MeV for both the 1999 
(open diamonds) and the 2000 (crosses) LXeGRIT electronics configuration. The 
gain in 2000 was about twice the gain in 1999. Right: energy resolution versus 
energy, showing the dependence expected from Poisson statistic corrected by 

a constant term. 



5 22^^(0.511 and 1.275 MeV), i37cs(0.662 MeV), 88y(0.898 and 1.836 MeV), 
6°Co(1.173 and 1.332 MeV), ^0K(1.465 MeV), Am-Be (4.43 MeV). 
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2 Gamma ray tracking 

The most general interaction sequence useful for Compton imaging is a first 
Compton scatter followed by full absorption of the scattered 7-ray in one or 
more interactions 0. In LXeGRIT, a large fraction of multi-site events has only 
two interactions (Compton scatter followed by photoabsorption) . The fraction 
of events with more than 3 interactions is at most 5% of the fully contained 
events for energies below 5 MeV and is not considered in the following. The two 
cases for 2- or 3- interactions are treated separately in Sec. 12.11 and Sec. 12.21 

Events with interactions other than Compton scattering and photoabsorption 
are not considered here for Compton imaging. The relative fraction of these 
events increases with energy and becomes dominant above ~5 MeV. The cross 
section for pair production turns on at 1.022 MeV and, in LXe, equals the 
Compton cross section at ~6 MeV. Relativistic electrons lose energy both by 
ionization and by radiation of secondary photons (Bremsstrahlung), which, if 
of sufficiently high energy, may be detected at a separate location. For electrons 
the radiative energy loss overcomes the one due to ionization above the critical 
energy = 610 MeV / (Z + 1.24) for liquids and sohds [15j. For Xe, Z = 54 
and Ec = 11 MeV. Experimentally, we observe considerable modifications to 
the energy spectrum for the 4.4 MeV Am-Be calibration source. 



2.1 2-site events 

In the case of 2-site events, assuming no prior knowledge the right sequence 
can be guessed with a 50% success rate. This improves using the energy shar- 
ing between the two interactions that is, in many cases, highly asymmetric. 
For LXeGRIT, the argument goes like this: for energies larger than ~2 MeV 
a 7-ray is more likely to be stopped in the fiducial volume with only two 
interactions if it does lose a large fraction of the initial energy in the first 
interaction. If the energy lost in the first interaction is small, the scattered 
photon will most likely interact more than once before being absorbed, i.e. 
would be classified as a 3(+)-site event. Fig. [7] shows the Ei and E2 energy 
spectra for 0.662 MeV (^^^Cs), 0.898 and 1.836 MeV (^^Y) photons. The ^^^Cs 
source was collimated to a beam with a lateral spread of ~3 mm, which al- 
lows us to tag the first interaction as the one within the collimator aperture. 
The ^^Y source was at a distance of 2 m above the detector, without any 
collimation. Since the source position is known, it is possible to use Compton 
imaging (Sec. [3]) and track each 7-ray assuming the two possible sequences. 
The sequence which gives the correct source position is then chosen as the true 



^ For additional reading on 7-ray tracking, see e.g. |ll|12|13|14j 
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one. For 1.836 MeV the Ei and E2 spectra almost mirror each other (Fig. [7]- 
right), and the situation Ei > E2 is clearly the most likely. At lower energies 
(0.662 and 0.898 MeV) the two spectra are much more similar and the E1/E2 
asymmetry is no more a good argument. The minimum in the Ei spectrum, 
clearly visible for all the three energies, corresponds to 90° scatter angle (the 
corresponding Ei = Ef^J {rrieC^ + Etot) is marked with a vertical dashed line). 
It is a geometrical artifact, due to the direction of the incident 7-rays along 
the detector z axis, such that for a 7-ray scattered at 90° Zi — Z2, while a min- 
imum separation of about 3 mm along the z-axis is required to ensure a good 
energy determination [2]. In Fig. ^left the ratio < Ei > /E^ is plotted vs. 
Ej for MC data, where < i?i > is the mean of the first energy deposition and 
E^ is the nominal energy of the 7-ray. The trend is quite clear: < Ei > / E^ 
increases with E^ and < ^1 > / E^ > 75% for E^ > 2 MeV. Fig. E-left is 




0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.5 1.0 1.5 2.0 

energy [MeV] energy [MeV] energy [MeV] 



Fig. 7. El, E2 distribution for fully contained 2-site events. From left to right: 
0.662 MeV (^^'^Cs), 0.898 and 1.836 MeV (^^Y). 




Fig. 8. Left: < Ei > / E^ vs. E^, from MC data. Superimposed f{E^) as defined in 
Eqs. [2]l3] Right: Efficiency and contamination in sequencing 2-site events according 
to the procedure described in Sec. 12.11 



suggestive of an optimized selection on the energy sharing, requiring the ratio 



8 



Ei/Etot to be larger than some value f[Etot) 

f(Etot) = 0.85 (^1 - ; Etot > 2.4 MeV (2) 

f{Etot) = 0.5 ; 1 MeV < Etot < 2.4 MeV (3) 

shown in Fig. ^left with Etot = E^. Fig. fright gives efficiency and con- 
tamination for this sequencing procedure in the energy range 1-10 MeV. The 
efficiency is as high as 86% at 2 MeV and saturates (>98%) at 5 MeV. Below 
1 MeV no optimum criterion for sequencing 2-site events has been found. 

2.2 Multi-site events 

When more than two interactions are available, the time sequence is, in prin- 
ciple, univocally determined by Compton kinematics. In the general case of a 
7-ray which undergoes — 1 Compton scatters and is photoabsorbed in the 
N^^ interaction energy and momentum conservation is written as 

EI, = EJ + E'r ■ pU=p1+P- (4) 



with E] {i = 0, . . . , N — 1) and E"? {i = 1, . . . , A^) the energy of the 7-ray 
and the scattered electron after interaction i; Eq is the energy of the incoming 
photon and Pi are the corresponding momenta. The electron scatter angle is 
not measured and is ignored in the following. For the photon scatter angle cp 



1 1 E 
— - , with: Wi = 

Wi Wi+i moc^ 



cosipi = l + — -— — , with: PV^ = — ^ (5) 



For a given interaction sequence, the interaction locations determine geomet- 
rically N — 2 photon scatter angles (fgeo i = 2, . . . , 

11 i ' U z-l-1 / \ 

COS(^geoi = |-> 11^ — r (6) 
I Ui\\Ui+i\ 



where iti = (xj - i/i - Zi - Zi_i). 

N — 1 Compton scatter angles are measured by the energy deposits accord- 
ing to equation [5l noting that Ej = Y.jLi+i Ej {i = 0, . . . , N — 1). This redun- 
dant information allows testing of the sequence of the interaction points based 
solely on kinematics. A straightforward test statistic consists in summing the 
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differences of the scatter angles quadratically, weighting the summands with 
the measurement errors: 



rp _ 1 (cos i^i - COS yPgeo if . . 

with: a] = a^^^ + a^^^ ^^^^ ^ 

Ideally, the test statistic would be zero for the correct sequence if the photon is 
fully contained. With measurement errors, is always greater than zero, but 
the correct interaction sequence is still most likely to produce the minimum 
value of the test statistic. A straightforward interpretation as a reduced 
distribution is not possible due to the non-Gaussian shape of the probability 
distribution in ip. For each triplet of interactions (Tcosip and (Xcosvpgco then are 
computed 



E 

fc=l 



Ui+l\ 



Uj^k cosy^gc 
I w'iP 



Mi+l,fc COSV^gco 



4 



\\Ui\-\Ui+i\ \Ui+l\ ) 

with: k spatial coordinate index and 

CTfc position uncertainty on each coordinate 

1 \' 



cos <p,i 



(9) 



We consider here only the case of 3 interactions, by far the most likely. There 
are six (3!) possible sequences to start with, i.e. assuming no additional knowl- 
edge the right sequence is chosen 17% of the times. The efficiency of this proce- 
dure is shown in Fig. O-Ze/t, using MC data. It is ~55% at 2 MeV and exceeds 
60% above 5 MeV. In LXeGRIT the energy resolution is the limiting factor. 
The fraction of wrongly sequenced events {contamination) is also shown in 
Fig. [9]-/e/t. In this case, since no event is rejected, it is just the complement 
to 1 of the efficiency. It is worthwhile to note that the most frequent confusion 
of the interaction sequence involves the swap of second and third interaction, 
while the first interaction is properly found. Such events, counted here under 
"confusion" , remain usable for imaging if the separation between second and 
third interaction is considerably shorter than the separation of first and second 
interaction. In this case, the wrong sequence leads to tails in angular resolu- 
tion but not to a "conversion" of source photons into background photons. 
Assuming the efficiency shown in Fig. [9]-/e/t as an upper limit, the algorithm 
performance can be improved applying further selections a posteriori, with 
the goal of keeping the efficiency as close as possible to the one in Fig. ^left 
while reducing the contamination. 
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The most powerful variable is Ei/Etot, where Ei is the energy deposited in the 
yirsi interaction 0. Fig. [TO]shows Ei/Etot for the 0.898 and 1.836 MeV line 
(experimental data), when all the events are considered and selecting only 
events with the time sequence correctly identified, i.e. events for which the 
source is correctly imaged. At 0.898 MeV there is no clear correlation between 
Ei/Etot and finding the right sequence, but at 1.836 MeV a large fraction 
of wrongly reconstructed events shows up at Ei/Etot < 0.3. The impact of 
selecting Ei/Etot > 0.3 has been studied over the energy range 0.5-10 MeV 
using MC data. The efficiency of the 7-tracking procedure combined with this a 
posteriori selection is shown in Fig. fright together with the contamination, 
here defined as the fraction of events (in the full energy peak) which are 
wrongly sequenced and have Ei/Etot > 0.3. This technique works well above 
2 MeV, while at lower energies the reduction in efficiency combined with poor 
rejection power makes it counterproductive, as also seen in Fig. HUlleft. 




Fig. 9. Efficiency and contamination in the reconstruction of the correct time se- 
quence for 3-site event, as calculated using MC data and the actual algorithm used 
for experimental data. Left: without applying any further selection. In this case 
events can only be correctly or wrongly sequenced, since no event is rejected. Right: 
selecting Ei/Etot > 0.3 a posteriori. The contamination fraction is here defined as 
the fraction of events wrongly sequenced which have Ei/Etot > 0.3. 



3 Compton imaging 



3.1 Angular resolution 



If the source position is known, two independent measurements of the first 
Compton scatter angle {(p and '^gco) are given. The difference between the two 

^ First according to the 7-tracking algorithm 
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Fig. 10. El/ Etot for 0.898 {left) and 1.836 {right) MeV 7-rays. From experimental 
data, imposing full energy containment. 

gives the measure of the angular resolution (ARM). In this section ARM and 
1-0" angular resolution are interchangeable. Three sources of error limit the an- 
gular resolution: energy resolution, which limits the precision in measuring 
position resolution, which limits the precision in measuring v^geo! and Doppler 
broadening. 



The uncertainty on the scatter angle due to the energy resolution is: 




Ei{Ei + 2E2)AE2 



ELEl 



(10) 



where Etot is the initial energy of the 7-ray, Ei the energy deposited in the 
first interaction, E2 = Etot — Ei; in Eq. [10] there are only two free parameters, 
e.g. (f and Etot- Its behavior as a function of (f and of Etot is shown in Fig. [TTl 
The curve for the (f dependence is obtained for Etot=^-836 MeV, the one for 
the Etot dependence is obtained integrating over (f < 60° according to the 
Klein-Nishina cross section. An energy threshold of 150 keV and an energy 
resolution of 10% / ,Je [MeV] (FWHM) have been assumed. 

The uncertainty on the direction of the scattered 7-ray, A(fgeo, assuming a 
separation between the first two interaction locations | t?| large compared to 
ax, o-y and az, is given by 

Ay^geo = |^7(A</Pgeo)^. + (A(^geo)^ + (A^^geo)^ (11) 




with (Av5geo)a = j^o-^^l - ] ,a = x,y,z. In Fig. [H] A^^geo is shown 

for =30 mm. 
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Fig. 11. Expected angular resolution for LXeGRIT; Aip and A(^geo have been 
obtained using Eqs. llOlllll and combined to give the final ARM according to Eg .[121 
Left: expected angular resolution vs. scatter angle if for a fixed energy (1.836 MeV). 
Right: expected angular resolution for LXeGRIT vs. energy, selecting ip <60°, i.e. 
forward scattering. 

The Compton formula in Eq. [1] gives the scattering angle if the incident pho- 
tons were to interact with stationary free electrons. Doppler broadening con- 
stitutes an irreducible limitation to angular resolution for a CT and its effect 
is larger for target materials of larger atomic number, such as Ge or Xe, com- 
pared to Si or liquid scintillators [16\. For low energy (few 100 keV) 7-rays, 
the uncertainty in the scatter angle due to Doppler broadening contributes 
significantly to the overall ARM. Once energy resolution and position resolu- 
tion are taken into account, Doppler broadening plays a rather negligible role 
for LXeGRIT (Fig. [TT]) . which is designed to image 7-rays of energy 0.5 MeV 
or higher. 

Neglecting Doppler broadening, the overall angular resolution A(p is defined 

as 



Alp 



A^l, + Acp^ 



(12) 



At 1.836 MeV the 1 a angular resolution is about 3° for scatter angles up to 
60°, improving for more forward scattering. The dependence of Aip on the in- 
teraction separation is shown in Fig. lT2lriaht. for different position resolution. 
Given the typical separation of the order of few cm shown in Fig. W2lleft. a 
mm position resolution is required for a good imaging performance. 

The ARM spectra for real data, 1.836 MeV 7-rays, is shown in Fig. [T31 For a 
realistic comparison with expectation, the analysis is based on MC data. In 
this way it is also possible to separate the response for 2- and 3-site events. 
The result for the energy band 0.5-10 MeV is shown in Fig. [TU-Ze/t. In Fig. UM 
right MC data and experimental data are compared for the lines: 0.662 (^^''Cs), 
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Fig. 12. Left: 3D separation between the first and the second interaction for 
1.836 MeV fuhy contained 7-rays. Right: Angular spread Aipg^o vs. 3D separation for 
different values of the position resolution. LXeGRIT achieves a position resolution 
of 0.85 mm. 



0.898 (8^Y), 1.275 (^^Na) and 1.836 
and 1.836 MeV (3-site events). 



*Y) MeV (2-site events), and 0.898, 1.275 




Fig. 13. ARM spectra for 1.836 MeV 7-rays (^'^Y source). The standard deviation 
has been obtained fitting the ARM spectra with a Gaussian function. Left: 3-site 
events. Right: 2-site events. The dashed line indicates events with (p restricted to 
less than 70°. This selection makes the ARM distribution sensibly narrower, getting 
rid of the extended tails due to large scatter angles. The standard deviation has 
been derived from this selected sample. 

The angular resolution vs. ip is shown in Fig. [13, from the same data (2- and 
3-site events combined); the experimental points are compared to the expected 
angular resolution (as shown in Fig. UTl-left). showing a good agreement. The 
expectation for a position resolution degraded to 2 mm is also shown; for small 
(f the overall performance is compromised. 
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Fig. 14. ARM spread vs. energy obtained using accurate MC data and experimental 
data. The 2- and 3-site samples are shown separately. Left: MC data; the point at 
0.5 MeV in the 3-site data and the one at 10 MeV in the 2-site data have been 
omitted because of very little statistical significance. Full energy containment has 
been imposed. Right: Experimental data: 0.662 (^^^Cs), 0.898 (^^Y), 1.275 (^^Na) 
and 1.836 (^^Y) MeV for 2-site events, 0.898, 1.275 and 1.836 MeV for 3-site events. 
The corresponding MC curves have been superimposed. 




20 40 60 80 100 120 140 
Scatter Angle cp [°] 



Fig. 15. Angular resolution vs. scatter angle for a sample of 1.836 MeV 7-rays, in ip 
bins of 20°. Superimposed, the expected angular resolution assuming two different 
values for the position resolution. The expected LXeGRIT position resolution of 
about 0.85 mm agrees well with the data. A 2 mm position resolution, while still 
in good agreement with the data for scatter angles larger than 60°, is clearly ruled 
out by the two points at 20° and 40°, the only ones actually sensitive to Atpgeo. 

3.2 Maximum Likelihood Lmage Reconstruction 



The imaging problem deals with the derivation of the intensity distribution of 
the object region from the observational data, which, for a Compton telescope, 
can be represented in a 3D binned data space consisting of the scatter direction 
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and the scatter angle y?: 

D^ = J2R^jfj + bi + Ni (13) 
j 

where Di is the observed number of counts, 6j the expected background, and 
Ni the statistical noise in the i^^ bin in data space, Rij is the instrumental 
response, i.e. the probability to detect a photon from the j^^ pixel in the 
object region in the i^^ bin in the data space, and fj is the flux in the j*'' 
pixel in the object region We restrict ourselves in the following to maximum 
likelihood fitting of single or few point sources (plus background), scanning 
and testing a grid of image pixels for one source at a time. The probability 
of the observed data under a specific model {fj, bi} is given by the likelihood 
function L, defined by multiplying the probability of each bin 

^=n^^=n^(Ai/„&.) (14) 



As a counting experiment with fixed observing time the statistics in each bin 
is given by the Poisson distribution 

Di 

Pi = ^ e-^' for uji>Q 

Pi = l for tUi = 0, A = 

Pi = for tUi = 0, A > (15) 

where a;, = Pijfj + h is the expected number of counts in the i*^ bin and 
Pi the probability of having Di counts in the i*'' bin, given Ui. 

Taking the logarithm of the likelihood function, one obtains 

logL = ^ Alogo;, + C (16) 



where C = — X]jlog(AO is a constant with respect to the parameters fj 
and is therefore model independent and can be neglected. Maximizing L with 
respect to the flux distribution {fj} for the intensity fj results in the following 
set of equations: 



This is the general expression for maximum likelihood in binned mode. A pos- 
sible choice to solve the system in Eq. [T7] is the Newton- Raphson algorithm 
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[T7] , which separately estimates each pixel on the source parameters. The sta- 
tistical significance is then obtained from —2 log A, where A is the likelihood 
ratio of the two hypotheses background only and background plus source [18j. 
—2 log A follows a xl distribution for an unknown point source and a Xi dis- 
tribution for a known source. 

A different approach, known as list mode likelihood method, aims at recon- 
structing an image on an event-by-event basis. This method may be derived 
from the binned likelihood method by increasing the number of bins in the 
data space until each bin contains either or 1 event. If only bins that contain 
an event are considered, Eq. [T7] deals with events rather than bins. The prob- 
ability to detect a photon from the j*'' pixel of the object region in the i^^ bin 
turns into the probability density for the i*^ photon in the data space. The 
list mode maximum likelihood method is very useful in overcoming problems 
of storage in computer memory and in reducing the CPU time needed for the 
calculation for cases of sparsely populated dataspaces, e.g., when each event 
contains multiple parameters relevant to the imaging problem. A detailed de- 
scription of this technique is given in [19]. 

In the case of a CT, the equation in list mode can be derived directly from 
the one in binned mode. One defines the source position {xo, ipo) in a 3D 
data space ((^, i^): where is some reference frame, e.g. longitude and 
latitude or right ascension and declination (Ra, Dec). The instrument response 
is equivalently given in a 3D data space as R^^\x^ i^i ^) or in a 2D data space 
as R^'^^if, if geo)- In the 3D data space R^^^x^i^^ 'flxo^^'o) has a conical shape 
with a half-opening angle of 45°and the vertex at the source location {xo, 
ipo) (Fig. [2]). _R(^)((^|99geo) is given by the probability of measuring Lp for a 
given scatter angle (pgeo', an example for LXeGRIT is shown in Fig. [111 The 
probability distribution is enhanced along the diagonal, i.e. for ip = v^gco, which 
is equivalent to having the ARM peak at 0°. R^^^ and i?*^^-* are connected 
through the relation 

Lp\Xo,i'o) 27ism(fgeod<fgeo = R^'^\^\^geo) cosipdxdip (18) 

We have presently implemented the list-mode likelihood method only without 
the background term, i.e. the expected number of counts in the i^^ bin is now 




(19) 



j 



The logarithm of the likelihood function is written as 



logL = ^ Aiog(E4 /.) - 



f^ + c 



(20) 
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If the bin size is reduced until each bin has at most one count, Di can only be 
or 1. The first term in Eq. [20] is now the sum over events rather than bins. 
Introducing an index ie to indicate events and id to indicate bins, Eq. [2D] can 
be rewritten as 

log L = ^ log - E R^jf^ +C (21) 

ie j id,j 

Approximating R^'^^Lp, v^geo) by a Gaussian G{lp, V^) for each value of <^ 

i?(^)(^,^geo)=/^G(^,V^^)A^^,„ (22) 



and replacing R^^^ in Eq. [5T] according to Eq. [TS] 

log L ^ E log f E ^ V E log ^^^"1^ ^ - E /. + C (23) 



where Tj = J2id ^fdj 1^ ^1^^ sensitivity to the j*'* pixel in the object region and 
the second term is a constant since the parameters J^, x? do not change for 
each event. Eq. [23] is therefore simplified as 

log = E log I E ^ I - E + C (24) 



which is the likelihood function in list mode. We are also applying a Newton- 
Raphson algorithm to maximize the list-mode likelihood function in equa- 
tion [M 



4 Imaging results 



The Maximum Likelihood imaging techniques in list mode has been used to 
produce images of calibration 7-ray sources. Fig. [17] shows the energy spec- 
trum obtained from exposing LXeGRIT to a 2738 kBq ^'^Y source at a distance 
of 2 m, on axis, for about 90 minutes; 3-site events have been selected. Be- 
fore any selection, the main features in the energy spectrum are the two 
lines (0.898 and 1.836 MeV), together with a continuum which extends up to 
~3.7 MeV due to partially absorbed 7-rays and, above 1.836 MeV, to pile-up 
of independent 7-rays. The energy spectrum after selecting events in the ARM 
peak (a selection also called software collimation) has been superimposed. The 
continuum is reduced by a factor of 4 at 1.5 MeV and to a negligible fraction 
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Fig. 16. R^'^\(p\ipgco) or point-spread function (PSF) for LXeGRIT, as obtained 
for MC data for a Crab-like source (i.e. with a power law energy spectrum with 
index 2 in the energy band from 1 to 10 MeV) 25° off-axis. The PSF is given by 
the probability of detecting (phi_bar) for each c/Pgeo (phi_geo) . Left: 2-site events. 
Right: 3-site events. 

above 2 MeV. The intensity of the 1.836 MeV hne is reduced by 45% by the 
ARM cut, consistent with the results presented in Sec. [31 The z and energy 
distributions for each of the three interactions are shown in Fig. [181 for events 
in the 1.836 MeV full energy peak and after software collimation. The same 
distributions from MC data reproducing the experimental conditions have 
been superimposed. The shape of the z distributions are as expected for a 
source on top of the detector, given that the first scatter is most likely in the 
forward direction. The image of the source for the 1.836 MeV line is shown 
in Fig. [ini reconstructed with a list mode Newton-Raphson algorithm. The 
source location is correctly determined with an accuracy of about one pixel, 
i.e. r. 

A second example is the resolved image of two calibration sources, ^°Co (1.17 
and 1.33 MeV) and ^^Na (1.27 MeV). The two sources were placed ~1.7 m 
above the detector with angular separation of ~10°. A fiat diffuse background 
and 100% detection efficiency Tj for each pixel in the object region were as- 
sumed, together with a variance of V^=3.5°. 



Conclusions 

LXeGRIT is the first fully developed and tested prototype of Compton Tele- 
scope based on a single position sensitive detector, such as a LXeTPC with 
combined event energy and 3D localization in one large homogeneous volume. 
In this paper the details of its performance in imaging MeV 7-ray sources have 
been presented. LXeGRIT has shown good performance as a 7-ray imager. 
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Fig. 17. Y 3-site energy spectrum before and after software collimation. 





el [MeV] e2 [MeV] eS [MeV] 

Fig. 18. Multi-site events selecting the 1.836 MeV line and applying software col- 
limation. Top: z distribution for the first, second and third interaction. Bottom: 
energy spectra for the first, second and third interaction. 



achieving an angular resolution of ~4° at 1.8 MeV, consistent with expecta- 
tions based on energy resolution, position resolution and geometry of its TPC 
Maximum Likelihood imaging techniques have been successfully applied to the 
LXeGRIT data. 
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Fig. 19. Maximum likelihood (list mode) image of an ^^Y source 2 m above the 
detector. 






Fig. 20. Maximum likelihood (list mode) resolved image of two calibration sources, 
^°Co and 22Na. 
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